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ABSTRACT 

Cas A is a Galactic supernova re mnant, whose supern ova explosion is observed to be of Type lib 
from spectroscopy of its light echo (jKrause et al.l 120081 ) . Having its SN type known, observational 
constraints on the mass loss history of Cas A's progenitor can provide crucial information on the final 
fate of massive stars. In this paper, we study X-ray characteristics o f shocked ambient gas of Cas A 
using the 1 Msec observation with the Chandra X-ray Observatory (|Hwang et al.ll2004l ). and try to 
constrain the mass loss history of the progenitor star. We identify thermal emission from the shocked 
ambient gas along the outer boundary of the remnant. Comparison of measured ra dial variat i ons o f 
spectroscopic parameters of the shocked ambient gas to the self-similar solutions of lChevalierl (|1982f ) 
show that Cas A is expanding into a circumstellar wind rather than into a uniform medium. We 
estimate a wind density nn ~ 0.9 ± 0.3 cm~ 3 at the current outer radius of the remnant (~ 3 pc), 
which we interpret as a dense slow wind from a red supergiant (RSG) star. Our results suggest that 
the progenitor star of Cas A had an initial mass around ~ 16 Mq, and its mass before the explosion 
was about 5 M , with uncertainties of a few tens of percent. Furthermore, the results suggest that, 
among the mass lost from the progenitor star (~11 Mq), a significant amount (more than 6 M ) 
could have been via its RSG wind. 
Subject headings: ISM: individual (Cas A) — supernova remnants — X-rays: ISM — stars: winds 



1. INTRODUCTION 

The mass of a star is the most crucial parameter that 
governs its structure, evolution and fate. Stellar masses 
vary throughout the life of stars: stars are subject to 
mass loss by winds and may experience mass exchange 
with their companion stars. For massive stars, in par- 
ticular, mass loss plays a cr itical role in their evolu- 
tion (e.g-. lMaeder et al1l2000fl . but this is poorly under- 
stood part of stellar evolution theory. Mass loss also 
profoundly impacts the eventual supernova (SN) explo- 
sion (Woo sley et al. 2002) and the consequent evolution 
of the supernova remnant (SNR). Massive stars, with 
their strong ionizing radiation and powerful winds, clear 
out a large area around them throughout their lifetime 
(e.g., a ra dius larger than 15 pc for stars earlier than B0; 
Chevalier 1999). This implies that the immediate envi- 
ronment of a massive star is thus determined by the mass 
loss history of the star itself. It further implies that when 
the massive star undergoes core-collapse and explodes as 
a supernova, the supernova shocks expand inside the cir- 
cumstellar structure. In other words, the structure and 
the evolution of "young" core-collapse SNRs reflects the 
pre-supernova mass-loss history of the progenitors, pro- 
viding observational constraints on the nature of their 
progenitor stars. 

Observations have established that stars with initial 
mass in the range 8-15 Mq undergo core collapse during 
their red supergiant (RSG) phases and produce Type II- 
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P SNe (|Smartt et alJl2009t ). On the other hand, the fate 
of more massive star s is still not well e stablished (e.g., 
see the discussion in iSmith et al.|[20TlT ). For example, 
stars with masses of 15 - 30 Mq may explode during 
the RSG phase, producing Type II-L, lib, or Iln SNe 
while stars with mass greater than 30 Mq may evolve 
into Wolf-Rayet stars and produce Type lb or Ic SNe. 
But observational constraints on these associations are 
indirect and weak. Also, the fate of massive stars may 
be complicated by binary interactions. 

For a star that undergoes a SN explosion during the 
RSG phase (or soon after it evolved past the RSG phase), 
its remnant will primarily interact with the surround- 
ing RSG wind. Existence of s uch an RSG wind has 
been suggested for SN 198 7A (Chevali er fc EmmerineJ 
[19891: ISueerman et~aT] [2005]) . The extent of the RSG 
wind would be determined by the location where the 
ram pressure of the wind equals the pressure of the sur- 
rounding medium. A radius of 5 pc is estimated for a 
canonical case (jChevalierl 120051 M = 5 X lO -5 M yr~ 1 , 
v w = lSkms" 1 , p/k = 10 4 cm~ 3 K). Thus, the in- 
teraction of the SNR with the RSG wind can last for 
thousands of years. The Galactic core-collapse SNR 
G292.0+1.8 is likely an example that is curre ntly inter- 
actin g with the RSG wind of its progenitor (jLee et al.l 
120101 LEE2010 hereafter). By analyzing the X-ray emis- 
sion from the outer shock region of this ~ 3000 year-old 
SNR with the Chandra X-ray observatory, we found that 
the remnant has been expanding inside a wind. The esti- 
mated wind density (njj = 0.1 ~ 0.3 cm -3 ) at the current 
outer radius (~ 7.7 pc) of the remnant is consistent with 
a slow wind from a RSG star. 

Cas A, one of the best studied Galactic SNRs, also 
appears to be interac ting with its RSG wind . At its 
age of about 330 yrs (Thorstcn sen et aLll200lL and ref- 
erences therein), its morphology and expansion rates 
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are consiste nt with a model SNR interacting with a 
RSG wind (IChevalier fc Oishil l200l . It is also found 
that the observed characteristics of the X-ray ejecta 
knots are consistent with Cas A expanding into a wind 
([Laming fc Hwang! 120031 : iHwang fe Laming! 12009t h On 
the other hand, the direct identification of shocked RSG 
wind has been limited. The remnant contains slow mov- 
ing shocked circumstellar clumps, called quasi-stationary 
flocculi (QSF). While earlier studies suggested that these 
could be fragments of the RSG shell swept-up by the 
later Wolf-Rayet wind, more recent studies favor that 
these are more likely de nse clumps in the RSG wind 
(|Chevalier fc Oishl 120031) . While observations have re- 
vealed detailed physical properties of these QSFs, the 
global characteristics of the RSG wind have not been 
well established so far. 

The characteristics of the wind surrounding Cas A, i.e., 
the mass loss history of the progenitor of Cas A, is par- 
ticularly pertinent as we know the supernova spectrum 
of Cas A at t he time of its explo s ion from the ligh t echo 
observations (jKrause et all 120081: iRest et al.ll2011[) . The 
SN spectrum of Cas A shows that it is remarkably similar 
to that of SN 1993J, a prototy pical Type lib SN w hose 
progenitor is likely a RSG star (jNomoto et al.lfl993h . 

In this paper, we study the circumstellar environment 
of Cas A and try to constrain the mass loss history 
of the progenitor star. The basic method is similar to 
the one we applied to G292.0+1.8 in LEE2010, although 
here we use the 1 Msec Chandra observation of Cas A 
(jHwang et al.ll2004l ) . The data are briefly summarized in 
Section^ In Section^ we describe our detection of ther- 
mal X-ray emission from the shocked ambient gas and 
study its spectral characteristics. In Section [4] we show 
that the observed characteristics of the shocked ambient 
gas are consistent with those expected from a shocked 
RSG wind. Then the implication of our results on the 
nature of the progenitor star (Section [5]) and cosmic ray 
acceleration (Section [5]) is discussed. Section [7] summa- 
rizes our results. 

2. DATA 

Cas A was the first light targe t of the Chandra X-ray 
Observatory (|Hughes et al.l l2b00), and has been observed 
multiple times with Chandra. The one million second ex- 
posure of Cas A was obtained with the ACIS-S3 d etec- 
tor on board Chandra in 2004 (jHwang et al.ll2004f ). We 
use this one million second Chandra observation in our 
study. The 1 Msec observation data consist of 9 sepa- 
rate ObsIDs (4634, 4635, 4636, 4637, 4638, 4639, 5196, 
5319 an 5320) with their exposure time ranging between 
^40 ks and ^170 ks. We reprocessed the level 1 event 
files to create new level 2 event files. We applied param- 
eters of the standard Chandra pipeline process, except 
that we turned off the pixel randomization. CIAO 4.2 
and CALDB 4.2.2 were used for all the reprocessing and 
analysis. 

3. ANALYSIS 

3.1. Imaging Analysis 

The X-ray emission from Cas A is dominated by emis- 
sion from shocked metal rich ejecta. Along the outer 
boundary of Cas A are seen nonthermal filaments, which 
are synchrotron emission f rom the TeV ele ctrons accel- 
erated at the shock front ((Reynolds 2008, for a recent 



review) . On the other hand, our study focuses on diffuse 
thermal emission from the shocked circumstellar gas that 
is expected to be fainter than the nonthermal emission. 
Although non-thermal X-ray filaments are seen along 
most of the outer boundaries, there are regions with weak 
non-thermal emission, especially in the northwest. The 
X-ray emission in these regions show no prominent limb- 
brightening, and we suspect that this is emission from 
shocked circumstellar wind. 

The identification of the X-ray emission from shocked 
wind in Cas A is complicated not only due to non-thermal 
filaments but also due to contamination from ejecta emis- 
sion. Furthermore, for the faint emission from shocked 
circumstellar gas, non-uniform backgrounds could be 
problematic. We consider two background components 
that are spatially varying, thus hindering the identifica- 
tion of the emission from the shocked circumstellar gas; 
they are out-of-time event [j along the CCD readout di- 
rections and a dust-scat tered X-ray halo around bright 
sources ()Overb cck 1965). In Cas A, both background 
components primarily originate from the bright ring of 
the shocked metal-rich ejecta, whose radius is about three 
quarters of the outer radius of the remnant. In the fol- 
lowing, we describe our attempt to remove these two 
non-uniform background components. We note, how- 
ever, that our intention of removing background com- 
ponents is to improve the "imaging" analysis of the faint 
emission from the outer shock region. For our spectral 
analysis, we use original data and the effect of these back- 
ground components will be minimized by using the local 
background spectra. 

To remove the contamination from out-of-time (OOT) 
events, we create OOT event files from level 1 event files 
following the method outlined in the M. Markevitch's 
script make_readout_bg3. In short, this method substi- 
tutes chip-y (the readout direction) values of each event 
in the original level 1 event list with random ones, and 
processes them identically to the original event files (we 
refer them as OOT event files). From these OOT event 
files, we create OOT images. Then, the OOT images 
are subtracted from the image created using the original 
event files scaled for the ratio of the readout time over the 
frame time (0.041/3.2), resulting in the OOT-subtracted 
images. 

The second background component considered is the 
dust-scattered X-ray halo emission, which is formed 
around an X-ray source du e to small angle scatter i ng by 
foreground dust particle s (|Mauche fc Gorensteinl 119891 : 
iPredehl fc Schmittl I1995D . Physical modeling of such 
emission, especially for an extended source, is compli- 
cated. It depends on the distribution of dust along the 
line-of-sight and the energy of the X-ray photons. In- 
stead, we adopt a simple approach, where the spatial 
distribution of the halo emission is modeled as Gaus- 
sian profile around the source. In other words, the 
dust-scattered halo emission is estimated by ordinary 
Gaussian-convolution of the source image. The width 
and the amplitude of the convolution kernel are fitted to 
result in a flat background. This understated approach 
gives a poor approximation of the real halo emission near 
the source. However, we find that this provides a rea- 

6 http : //exc . harvard . edu/proposer/P0G/html/chap6 . html#sec : trailed- im 
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Fig. 1. — (a) Unprocessed Chandra image (3.34-7 keV) of Cas A. The image is smoothed with Gaussian beam of FWHM=4.5"for better 
visualization, (b) Same Chandra image as in the left but with the out-of-time event subtracted (see the text for details) . (c) Same Chandra 
image as in the center but the estimated dust scattered halo emission subtracted (see the text for details). 



sonable approximation for the outer regions that we are 
interested in. 

The steps of background subtraction are illustrated 
in Figure [T] The dominant background contribution in 
the original image (Figure QJa)) is the OOT contribu- 
tion seen along the direction of northwest to southeast 
which is the readout direction. The OOT contribution is 
effectively removed in the OOT-subtracted image (Fig- 
ure Eh)) while the remaining halo emission is clearly 
visible. And we see that the halo contribution is signifi- 
cantly reduced in our final image (Figure [He)). We note 
that the results depend on the chosen energy band. We 
find that the energy range of 3.34 - 7 keV, which does 
not include any prominent metal lines, provides the best 
result. In softer energy ranges, the strong line emission 
from the shocked metal-rich ejecta makes the background 
subtraction less optimal. 

In Figure [21 we reproduce the Figure QJc) in false-color 
with different contrast. It clearly shows the existence 
of faint diffuse emission around the outer boundary of 
Cas A. The nonthermal emission is relatively weak in 
the western boundary where the emission shows a faint 
plateau. In contrast to the non-thermal filaments, lit- 
tle hint of limb-brightening is seen in this region. The 
emission without the limb-brightening is indeed what one 
would expect from the thermal emission of the shocked 




Fig. 2. — Background subtracted Chandra image (3.34-7 keV) 
of Cas A smoothed with Gaussian beam of FWHM=4.5". White 
lines are the areas where the radial profiles in Figure|3]are extracted 
and labels mark their position angles. 



ambient gas when the remnant is expanding into a wind, 
as demonstrated in the case of G292.0+L8 (LEE2010), 
and they are promising candidate emission features from 
shocked circumstellar wind. We note that this emission, 
in projection, often extends beyond the nonthermal fila- 
ments which marks the location of the outer shock front. 
To test if the faint emission beyond the nonthermal fil- 
aments is real and not some artifacts of our background 
subtraction, we extract radial intensity profiles in sev- 
eral regions from the original image and also from the 
background-subtracted image. White lines in Figure [5] 
are the regions where the radial profiles are extracted 
(these regions are associated with the regions where spec- 
tra are extracted in Section f3.2[) . and their radial profiles 
are shown in Figure [3] For each profile in Figure G2 we 
mark the locations of the outer shock front that we vi- 
sually identified from the background-subtracted image 
(Figure [2]) . We find that the profiles extracted from the 
original image do show changes in their slope across the 
proposed boundaries, indicating they are real. As will 
be discussed in the next section, spectra of these regions 
confirm their ambient origin. 

3.2. Spectral Analysis 

We extract spectra from several regions around the 
outermost boundary believed to be dominated by ther- 
mal emission of shocked ambient gas. The extraction re- 
gions are shown in Figure [4] and the corresponding spec- 
tra are shown in Figure [5] The spectra are extracted 
from the original event files, thus the processes described 
in Section 13.11 do not apply for the spectral analysis. 
Instead, the background spectra are extracted from the 
regions right outside the boundary, therefore the contri- 
bution of the OOT event and the halo emission to the 
source spectra is should be negligible. We selected re- 
gions with little contamination from the ejecta emission 
and the nonthermal emission. The spectra of these re- 
gions clearly show Silicon and Sulfur lines. The overall 
spectral slope in the hard part of spectra is also softer 
than that of the nonthermal filament (an example spec- 
trum of nonthermal filament is also shown in Figure [5] for 
a comparison). Therefore, we conclude that the X-ray 
emission we have extracted is dominated by a thermal 
component. 

The extracted spectra are fit with a single nonequi- 
librium ioniz ation p lane-parallel shock model (vpshock 
in Xspec v!2, Borkowski ct al. 2001) with varying abun- 
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Fig. 3. — Radial intensity profiles (E=3.34 ~ 7 keV) extracted 
along the several regions shown in Figure [2] The black lines are 
from original image and the gray lines are from background sub- 
tracted image. The left panels show the radial profile in a log-log 
scale. The right panels show the same radial profiles in a linear- 
linear scale for the region near the outer boundary (marked as gray 
boxes on the left panes). The locations of the outer boundary vi- 
sually identified from Figure[2]are marked as vertical dashed lines. 

dances. The spec tra are not binned, ins tead we use Chu- 
razov weighting (Chur azov et al.l 119961 ). The net pho- 
ton count for each spectrum is between 5,000 and 20,000 
photons. The best fit models are overlaid in Figure [5j 
We find that a single vpshock model adequately fits the 
observed spectra. The fitted abundances are similar to 
or less than the solar values. A similar low abundance 
pattern has been also found from the outer shock of 
G292.0+1.8 (LEE2010). Therefore we suggest that the 
faint emission we see is thermal emission from shocked 
ambient gas. In Table[TJ we summarize the fit parameters 
for each regions, sorted by their position angle (P. A.). 
Througho ut this paper, we ad opt the explosion center 
from Thorstcn sen et aLl (|2001l ) as the remnant center. 
It is possible that there is some contribution from the 
nonthermal emission, and any residual nonthermal com- 
ponent would affect the fit parameters. We note that 
fitting the spectra with additional powerlaw component 
(with its photon index fixed at 2.5 from the fit to nearby 
nonthermal filaments) decreases the fitted temperature 
while increasing the ionization time scale, but other pa- 
rameters remain consistent within their error ranges. 

In FigurelHl we compare spectral parameters of each re- 
gion as a function of their position angles. The absorbing 
hydrogen column densities and electron temperatures are 
given from the spectral fit. The densities are estimated 
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Fig. 4. — Regions where the azimuthal spectra in Figure \5\ are 
extracted. The regions are labeled by their position angles from 
the remnant center. 
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Fig. 5. — Extracted Chandra spectra of the outer faint regions 
in Figure [4] and their best fit models. The bottom spectrum is 
from the nonthermal filament also marked in Figure [4] shown for 
comparison. 

from the fitted volume emission measures. The volume 
corresponding to each region is calculated by pixel-by- 
pixel integration of the path-lengths over the projected 
area. In this calculation, we assumed a sphere centered 
at the a ssumed center of the remnant at a distance of 
3.4 kpc (jReed et al.lll995l) . The radius of the sphere is 
separately determined for each region as the projected 
distance of the farthest (from the assumed center of the 
remnant) pixel in the region. The absorbing hydrogen 
column densities seem relatively low in the north com- 
pared with ot her directions, which seems consistent with 
the results of iHwang fc Lamina (|2012[ ). On the other 
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TABLE 1 

X-RAY SPECTRAL PARAMETERS OF THE REGIONS USED IN THE AZIMUTHAL ANALYSIS. 
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t Position angles to each region are measured east of north from the geometrical center. 
* For all individual spectra, the degree of freedom (v) is 454. 
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Fig. 6. — Azimuthal distribution of the absorbing column densi- 
ties (top), electron temperatures (middle) and the postshock den- 
sities (bottom) from the spectral fits of the outermost regions in 
Figure [5] 

hand, we could not see a clear azimuthal trend in the 
electron temperatures and the hydrogen densities. Sim- 
ilarly, Figure [7] shows variation of the electron tempera- 
tures and the hydrogen densities of each region as a func- 
tion of their distance from the remnant center. No clear 
indication of radial trends is noticed although the radius 
range that we are sampling is rather limited (2.7' < R 
< 3.2')- 

As has been demonstrated in the case of G292.0+1.8 
(LEE2010), one can investigate the radial structure of 
the "shocked" ambient gas to test whether the gas is 
of interstellar origin or of circumstellar origin. We note 
that the radial structure we are considering here is qual- 
itatively different from what is shown in Figure [3 which 



Fig. 7. — Electron temperatures (top) and the postshock den- 
sities (bottom) from the spectral fits of the outermost regions in 
FigurePolare plotted as a function of their radial distances from the 
center. If Cas A is expanding inside a wind, the radial variation 
of postshock density in this figure will follow the preshock density 
profile scaled by the shock compression ratio. The solid line is 
the best fit model with a simple wind profile (p oc r -2 ) and the 
gray area represents the formal error of the fit for 90% confidence 
interval. 



represents characteristics of "multiple" shocks of differ- 
ent directions (and at different radii from the center) . By 
radial structure, we mean radial variations of gas charac- 
teristics behind "a" shock, which represents accumulated 
history of shock propagation. Investigating the nature of 
ambient gas from the radial structure behind a shock re- 
quires a radial series of spectra that covers large enough 
radial range to test different models. We select regions 
in the northwestern boundary where emission from the 
nonthermal filaments is insignificant. Figure [8] shows the 
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regions where we extract the radial series of spectra. We 
allow the area of the regions to increase outwardly to 
compensate for decreasing surface brightness. Again, the 
background spectra are extracted from the regions right 
outside the shock boundary. In the inner region away 
from the shock boundary, the contribution of the OOT 
event and the halo emission may remain in the source 
spectra. However, we do not think that significantly af- 
fects our analysis because, as will be shown later, the 
emission from the CSM itself is much strong in this in- 
ner region. The spectra of individual regions are first fit 
independently as in the azimuthal analysis. The fitted 
metal abundances are found to be less than those of so- 
lar values, even in the bright innermost region well inside 
the remnant boundary. This suggests that most emission 
(if not all) of our selected regions is from shocked ambient 
gas. The fitted absorbing column densities (-/Vh) of these 
regions were consistent within their uncertainties. We 
consider that the variation of ./Vh in this small patch of 
the sky is negligible compared to uncertainties in the fit 
parameters and we refit the whole set of spectra with Ah 
fixed at the value of 1.2 x 10 22 cm -2 , an average from the 
independent fits. Figure [9] shows a selected sample of the 
spectra from these regions and their best fit models. And 
the fit results are summarized in Table [5] and plotted in 
Figure [TUJ It is clear that the observed emission measure 
monotonically increases toward the remnant center. On 
the other hand, the variation of the electron temperature 
is complicated. In the inner regions (radius smaller than 
2.5'), we see clear tendency of inwardly decreasing elec- 
tron temperature. But the temperature remains more or 
less flat in the outer regions (radius between 2.5' ~ 3.0'). 
For some elements (e.g., Ne, Mg, and Si), their abun- 
dances are higher for the inner regions. However, it is 
not clear whether this trend is real or is due to some sys- 
tematic uncertainties. The radial structure behind the 
shock will be further discussed in the next section. 
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Fig. 8. — Zoomed view of the northwestern part of Figurcl2lwith 
regions for the radial spectral analysis are marked as white boxes. 



4. NATURE OF THE AMBIENT GAS 

The radial structure of the shocked ambient gas can 
trace the radial density structu re of the medium into 
which the remnant is expanding. IChevalierl ([1982) found 




Energy (keV) 

Fig. 9. — Sample of Chandra spectra extracted from regions 
for the radial analysis in Figure [8] The uppermost spectra corre- 
sponds to the innermost region, and the lowermost spectra to the 
outermost region. Their best fit models are also shown. 

self-similar solutions for the interaction of expanding stel- 
lar ejecta (p ex r~ n , n > 5) with an external medium 
(p oc r~ s ) where s and n are the power law index of 
the density profile for ambient gas and for the ejecta, re- 
spectively: a remnant expanding in a uniform ambient 
medium (s = 0) would result in an inward-decreasing 
density and an inward-increasing temperature profile of 
shocked ambient gas, while a remnant expanding in a 
medium with a wind profile (s = 2) would show an 
inward-increasing density and an inward-decreasing tem- 
perature. The primary effect of n (the density structure 
of the ejecta) to the solution is change in the distance be- 
tween the forward shock and the contact discontinuity, 
and does not fundamentally change the overall behavior. 
Therefore, by studying the current radial profile of the 
density and temperature of shocked ambient gas, one can 
infer whether the SNR has been expanding in a density 
gradient. The utility of this method has been success- 
fully demonstrated in our study of Galactic core-collapse 
SNR G292. 0+1.8 (LEE2010). 

In Figure I10( we overlay model predicted radial struc- 
ture of the temperature and emission measure on top 
of the observed ones. We adopt self-similar solutions of 
Chevalier (1982) and project them assuming a spherical 
geometry. As in LEE2010, the projected temperature is 
taken as an average temperature along the line of sight 
weighted by density squared, ignoring the variation in the 
intrinsic emissivity of the gas at different temperatures. 

The increasing emission measure toward the remnant 
center, as can be seen from Figure [TOl is the behavior ex- 
pected for the remnant expanding inside a wind and we 
find that the observation is consistent with self-similar 
solutions with s = 2 and 6 < n < 7. The temperature 
variation is not consistent with any of the models, al- 
though the decrease in temperature in the inner part (R 
< 2.5 pc) is what is expected in the s = 2 models. As will 
be discussed below, direct comparison between observed 
temperature to model ones is difficult due to complicated 
shock physics, especially near the shock front. On the 
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TABLE 2 

X-RAY SPECTRAL PARAMETERS OF THE REGIONS USED IN THE RADIAL ANALYSIS. 
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Fig. 10. — Fitted electron temperatures and emission measures 
of the radial series spectra of Figure [9] plotted as a function of 
their radial distances from the center. Lines are model predictions 
from Chevalier (1982; see Section [4] for more details) for different 
sets of parameters, where s and n are power-law index for the 
density profile of the ambient medium and of the expanding ejecta, 
respectively. The temperatures and emission measures of the model 
predictions have been scaled to match the observed values. The 
inset of the bottom panel shows same data with different scales. For 
the models of s = 0, contact discontinuity is marked by cusps in the 
temperature profile. For the models of s = 2, contact discontinuity 
is where the temperature drops to zero and the emission measure 
diverge to infinity. 

other hand, the fact that the observed variation of the 
emission measure is in good agreement with models of 
s = 2 clearly suggest that the density of the ambient 



gas radially decreases, which confirms that Cas A has 
been expanding inside the wind. While our data suggest 
a power-law index n of ejecta envelop between 6 and 7, 
this may not be a reliable estimate as has been discussed 
in LE E2010. We note, however, that lLaming fe Hwanel 
(2003) studied ejecta emission in Cas A and found n be- 
tween 6 to 7. 

The temperatures from our X-ray spectral fits are elec- 
tron temperatures, while the temperatures of the models 
are ion temperatures. These two can differ considerably 
in non-radiative shocks faster than 1000 kms -1 (e.g., 
IGhavamian et all [2007, and references therein). The 
shock velocity measured from the proper motion of non- 
thermal filaments is about 5000 kms -1 (|DeLanev et al.l 
I2004D . while the measured electron temperature ranges 
between 2 to 2.5 keV (Table [1}. This corresponds to 
T e /Tp = 0.05, assuming all the shock energy has con- 
verted into the energy of the ordinary gas (i.e., no cosmic 
ray acceleration; however, see Section |5]). Previous stud- 
ies suggested that, for fast shocks (v s > 1000 kms -1 ), 
the temperature equilibration between electrons and pro- 
tons are low (T fi /T n < .1) and possibly scales with v~ 2 
(jGhavamian et al.l I2007J ). The observed electron tem- 
perature in Cas A and its inferred T c /T p , assuming no 
cosmic ray acceleration, seems to be consistent with ob- 
served values of other SNRs. Since it has been known 
that there is CR acceleration in Cas A, it is likely that 
the observed low electron temperature is partially due to 
an efficient cosmic ray acce leration at the forward shock 
(e.g., iVink fe Lamin g 2003) . Nevertheless, the measured 
electron temperature is likely not to trace the tempera- 
ture e xpected from the self-similar solutions of lChevalierl 
(1982) and the discrepancy between the models and the 
observed values likely represent incomplete physics of the 
self-similar models we used. However, while efficient cos- 
mic ray acceleration may affect some of our results, e.g, 
estimated preshock densities, we do not consider the ef- 
fects to be highly significant in the sense that they will 
not invalidate our primary result that Cas A has been 
expanding inside its wind. 

The properties of shocked ambient gas we estimated in 
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Section [3.21 are therefore the properties of wind material. 
The measured postshock densities in Figure [7] are from 
regions around the outermost boundary for di fferent po- 
sition angles, and the self-similar solutions of iChevalierl 
(1982) show that the density remains flat near the shock 
front. Thus, we consider these measured densities closely 
approximate wind densities at the location of the forward 
shock, scaled by the assumed postshock compression ra- 
tio. Fitting the observed densities in Fig. [7fb) with a 
wind profile assuming a constant postshock compression 
ratio of 4 gives a wind density profile of 



n H = 0. 



3pc 



The formal error from the fit (for 90% confidence inter- 
val) is about 5%. However, given the scatter of den- 
sity values for individual regions, which may reflect the 
asymmetric nature of the wind, an error of 30% would 
be more realistic. The inferred density profile of Cas A 
corresponds to the mass loss rate of 



M = 2.5 x 10~ 5 ( 



lOkms 



-)M Q yr- 1 . 



where v w is the wind velocity. The mass loss rate is too 
large to be that of the WR wind (M > 10~ 3 M Q yr" 1 for 
v w > 1000 km s~ ) and consistent wit h that of the RSG 
wind estimated from SNe studies (e.g.. lSramek fc Weilerl 
|2003|) . Therefore, our results support that Cas A is cur- 
rently expanding within the RSG win d. Using the mor- 
pholog y and expansion rates of Cas A, iChevalier fc OishU 
(|2003f) estimated explosion energy (E) and the total mass 
of the ejecta (M e ), but as a function of the wind den- 
sity. Substituting our observed wind density, we get 
E = 5 x 10 51 erg and M e = 4M Q . And the observed 
mass of the swept-up wind is around 6 M Q . For all these 
estimates, the uncertainties of our density estimate will 
propagate linearly, thus we expect about 30% of uncer- 
tainties. Note that these estimates may be further sub- 
ject to systematic uncertainties. For example, these val- 
ues are calculated assuming the RSG wind had a con- 
stant mass loss rate, and will be affected by a possible 
time variation of the mass loss rate. It was proposed 
that the dense wind of Cas A may not have extended di- 
rectly to the surface of the progenitor star and Cas A 
might have explod ed into a l ow density b ubble (e.g., 
Krause et a l. 2008) . Ho wever, iSchure et all (2008) and 
Hwang fc Lami ng (2009) found that the size of the bub- 
ble, if it existed, would not have been larger than 0.3 
pc, which is significantly smaller than the current radius 
of Cas A (~ 3 pc). Therefore, the effect of the bubble 
on ou r esti mates does not see m significant. IVink et all 
(|1996() and lFavata et all (|1997fl derived a total swept up 
mass ranging from 5 to 9 M© using the global X-ray 
spectrum of Cas A. While this is in agreement with ours, 
the mass estimates based on the fits to the global X-ray 
spectrum may suffer from uncertainties due to the spa- 
tial variation of X-ray spectra and the contribution from 
shocked ejecta. 

5. PROGENITOR OF Cas A 

A lower limit on the progenitor's initial mass can be 
estimated by combining all the masses associated with 
Cas A. If we combine the shocked wind mass (> 6 M Q ), 



the ejecta mass (4 M Q ), and the mass of the neutron star 
(~ 1.4 M©), we estimate the initial mass of the progeni- 
tor at least larger than 11 M Q . Another estimate of the 
initial mass comes from the SN type of Cas A. The SN 
spectra of Cas A observed by the light echos show that 
they can be classified as a Type lib SN (jKrause et al.1 
2008; Res t et al.l 1201 lh , which originate from the core 
collapse of massive stars that have lost most of their hy- 
drogen envelopes and consist of a nearly bare helium core 
at the time of its collapse. Thus the He core mass of 
the progenitor is supposed to be comparable to the stel- 
lar mass at the time of explosion, which is estimated as 
~ 5 M Q by summing the masses of the ejecta and the 
neutron star. And the He core mass of ~ 5 M© impl ies an 
initial mass of ~ 16 M Q (e.g.. lWooslev et al M2002f) with 
possibly a few tens percent of uncertainties. We note 
that this mass is close to the upp er limit of Type II-P 
SNe (~ 15M , iSmartt et al.ll2009l ). For a comparison, 
the progenitor star of the SN 19 93J, the prototypica l SN 
Type lib, was a red supergiant (jNomoto et al.lll993l and 
references therein) with a main sequence mass probably 
in between 13 and 20 M ? ([Wooslev et al.lH99l and ref- 
erences therein). 

Our analysis of the X-ray data suggests that Cas A 
has been expanding in a dense RSG wind whose swept- 
up mass is around 6 M . This large mass of swept- 
up RSG wind seems to indicate that the mass loss 
of the Cas A's progenitor might have been primarily 
via the RSG wind. However, stellar evolution mod- 
els, with c ommonly adopted m ass-loss rates for RSG 
stars (e.g., Ide Jager et al.l Il988t) . predict significantly 
less mass loss via the RSG wi nd for a single ^16 M© 
star (e.g., iWooslev et al.l [2002]) . For the progenitor of 
SN 1993 J, the mass loss is attributed to Roche lobe 
overflow in the binary system (jStancliffe et al.ll2009l and 
references therein). This was supp orted by the possi- 
ble existence of a binary companion (M aund et al.ll2004l : 
iMaund fc Smarttl 12009). A similar binary system als o 
has been suggested for Cas A (e.g.. lYoung et al.ll2006|> . 
although no companion star has been currently identi- 
fied. We note, however, that the mass loss rate of cool 
stars such as RSG stars is largely uncertain, and mod- 
els may have underestimated the mass loss during the 
RSG phase. For example, stellar pulsations may trigger 
an unusually strong wind which sub stantially in c reases 
the mass loss rate in the RSG phase. lYoon et al.1 (2010) 
suggested that, if that happens, a single 20 M Q star can 
explode as a Type lib SN after losing most of its hydro- 
gen envelope. 

6. COSMIC RAY ACCELERATION 

So far we have neglected effects of cosmic ray accel- 
eration in our analysis. However, the existence of X-ray 
nonthermal filaments around Cas A clearly indicates that 
electrons are accelerated up to the TeV energies. Obser- 
vations also suggest that magnetic fields are amplified at 
the shock, w hich suggests th e acceleration of cosmic ray 
protons (see [Reynolds 2008, for a review) . And modeling 
of the broadband nonthermal emission from the radio to 
the gamma-r ay bands indicates a cceleration of protons 
as well (e.g.. iBerezhko et al.ll2003[ ). Inclusion of cosmic 
ray acceleration will change the self-similar solutions we 
used, where a principal effect of efficient cosmic ray accel- 
eration would be a reduction of the distances between the 
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forward shock and t he reverse sho c k, as demonstrated in 
the Tycho SNR by iWarren et al.l (J2005D . However, the 
tendency of increasing density and decreasing temper- 
ature t oward the contact discon tinuity is not likely af- 
fected (jDecourchelle et all 12000). Thus, our conclusion 
that Cas A is expanding inside a wind would be valid 
even if cosmic acceleration were efficient in Cas A. On 
the other hand, we might have overestimated the wind 
density as the cosmic ray acceleration can increase the 
postshock gas compression ratio. For a strong relativistic 
shock, the compression ratio would be 7. The compres- 
sion ratio will be further increased if accelerated cosmic 
rays efficiently escape the system. 

We have assumed that the thermal emission we ana- 
lyzed is from the postshock gas. In particular, we as- 
sumed that those regions we used for azimuthal analysis 
represents postshock gas right behind the forward shock. 
However, most of these regions are located beyond the 
nonthermal filaments that are the locations of forward 
shocks. We note that the intensity of the nonthermal fil- 
aments varies along the boundary, and it is particularly 
faint in the western boundary where the thermal emis- 
sion is clearly seen. In this regard, we consider that the 
thermal emission beyond the nonthermal filaments are 
likely a projection of different shock fronts. It indicates 
that not all outer shocks of Cas A accelerate cosmic ray 
electrons efficiently enough to show X-ray nonthermal fil- 
aments. This further suggests that the mean expansion 
velocity of the shock front with nonthermal X-ray emis- 
sion is slower than those without nonthermal X-ray emis- 
sion. The slowing down of the shock might have been a 
consequence of efficient cosmic ray acceleration. In other 
words, there are regions where the acceleration has been 
less efficient, resulting in a larger shock radius that is 
sufficiently large to extend beyond the non-thermal fila- 
ments even in projection. This may be the case if the effi- 
ciency of CR acceleration depends on the angle between 
the shock no rmal and the magne tic field, as suggested 
for SN 1006 (jRevnosoet all 120131) . On the other hand, 
it might be that the preshock condition necessary for ef- 
ficient cosmic ray acceleration is related to the slowing 
down of the shock (e.g., density). Nevertheless, under- 
standing differences of these shocks may provide valuable 
insight on the cosmic ray acceleration in the shocks. 

7. SUMMARY 

Using 1 Msec Chandra observation, we present our de- 
tailed analysis of the blast wave in the Galactic SNR 
Cas A. By analyzing the X-ray emission from the shocked 



ambient gas, we have studied the nature of the ambient 
medium around Cas A, and also the nature of the pro- 
genitor star. Our results provide observational evidence 
of Cas A interacting with the RSG winds from the pro- 
genitor star. 

• After removing spatially- varying background com- 
ponents, wc identified faint emission, with no 
prominent limb-brightening, around the outer 
boundary of Cas A. 

• X-ray spectra of these regions are well fit with the 
emission from a thermal plasma of subsolar abun- 
dances, indicating they are shocked ambient gas. 

• In the western boundary where nonthermal emis- 
sion is relatively weak, we extracted a radial series 
of spectra. Comparison with self-similar solutions 
of IChevalien (1982) indicates that Cas A is cur- 
rently expanding inside a wind, rather than a uni- 
form medium. 

• The postshock density estimates from the spectra 
of the outermost shocked wind range between 3 and 
5 cm -3 , implying a preshock wind density of ~ 1 
cm~ 3 at the current outer radius of Cas A (~ 3 
pc). The high density of the wind suggests that it 
is a dense wind from a red supergiant (RSG) star. 



• Combined with results of lChevalier fc Oi shi (2003), 
our measurements suggest that the ejected mass 
during the SN explosion was about ~ 4M Q , im- 
plying total mass of the star (including the mass 
of the neutron star) at the e xplosion of ~ 5 M& . 
Since Cas A was Type lib SN (jKrause et al.l l2008). 
we estimate ~ 16 M© as the initial mass of the 
progenitor star. 

• Estimated explosion energy of 5 x 10 51 erg is larger 
than that of normal core-collapse SNe (~ 10 51 
erg), although not u nprecedented for Type lib SNe 
(Ha muy et al.l 120091 ) . On the other hand, the ex- 
plosion energy would have been overestimated if 
wind density is overestimated. 

• Our results suggest that the swept-up mass of the 
RSG wind is about 6 M Q . This indicates that the 
progenitor of Cas A may have lost its mass primar- 
ily via the RSG wind. 
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